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Abstract. Global Positioning System (GPS) satellites 
broadcast at frequencies of 1,575.42 MHz (LI) and of 
1,227.60 MHz (L2). The dispersive property of the iono- 
sphere can be used to combine independent measurements 
at the two frequencies to estimate the total electron con- 
tent (TEC) between a GPS receiver site and a broadcast- 
ing satellite. Such measurements, made at sites near to 
Very Long Baseline Interferometry (VLBI) sites, can be 
used to estimate the ionospheric contribution to VLBI 
observables. For our 1991.9 astrometric VLBI experiment 
in which we obtained group-delay observations in the 8.4 
and 2.3 GHz bands simultaneously, we found that the GPS 
and VLBI determinations of the ionosphere delays agreed 
with root-mean-square differences below 0.15 ns for inter- 
continental baselines and 0.10 ns for continental ones. We 
also successfully applied the GPS-based procedure to re- 
duce the ionospheric effect in phase delays used for high 
precision differenced astromctry at 8.4 GHz for this same 
experiment. 
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1. Introduction 

VLBI provides unprecedentedly accurate angular resolu- 
tion through observations of celestial bodies with radio 
telescopes spread over the Earth's surface. Each observ- 
ing station records data on magnetic tapes. The local- 
oscillator signals and the time-tagging of the data are gov- 
erned by hydrogen-maser frequency standards. The tapes 
are processed in special-purpose correlators to determine 
the so-called VLBI observables: group and phase delays 
and phase-delay rates. 

A main problem in determining the sky positions of 
celestial radio sources from these VLBI observables is the 
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effect of the Earth's ionosphere on them. The use of GPS 
satellite data to remove this effect forms the thrust of this 
paper. 

The ionosphere is characterized by its content of free 
electrons and ions. The F2 layer of the ionosphere has 
the largest density of charged particles, with values up to 
3 X 10^^ m~^. The total electron content per square meter 
along a line of sight is the number of electrons in a column 
of one square meter cross section along the ray path: 

TEC = / N -dh, (1) 
Jo 

where N is the spatial density of electrons, h is the co- 
ordinate of propagation of the wave, and ho corresponds 
to the effective end of the ionosphere. TEC is highly vari- 
able and depends on several factors, such as local time, 
geographical location, season, and solar activity. TEC can 
have values between 1 TECU (or TEC unit, defined as 
10^^ m~^) and 10^ TECU. Epochs of greater solar activ- 
ity cause higher values of the TEC. 

The ionosphere affects the phase and group delays op- 
positely (to first order, see, e.g., Thompson et al. 1986| ): 

Ar = T^-TEC, (2) 



where 



40.3 m'^s 



is the speed of light (m-s ^), 



and v the frequency (Hz), and where we neglect magnetic 
field effects and assume ly is large compared with the local 
plasma frequency (for an extreme case, the plasma fre- 
quency is of 15 MHz) . The negative sign applies for phase 
delays and the positive sign for group delays. In standard 
astrometric VLBI experiments observations are made si- 
multaneously in two well separated bands of frequencies 
in order to estimate the ionospheric effect. A nearly vac- 
uum equivalent delay can be obtained from the following 
expression: 



(l^l/l^2)^Tl - T2 



(3) 



(^^iM)2-i ' 

where Tj is the delay -either group or phase- at frequency 
i^i {i = 1,2, i^i > V2). 
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Thus, with dual frequency observations, the ionosphere 
effect can largely be removed from the VLBI data. Such 
removal can also be made for single-frequency observa- 
tions via Eq. (|2|), if estimates of the TEC along the lines 
of sight of the radio telescopes are available from other 
observations. 



Guirado et al. (1995) showed that it is possible to esti- 
mate the ionospheric effect with accuracy useful for astro- 
metric purposes from Faraday-rotation measurements. In 
this work, the authors used a "clipped" sinusoidal func- 
tion to model the diurnal behavior of the TEC. In this 
model, the night component is constant and equal to the 
minimum TEC value, and the day component is expressed 
as the positive part of a sinusoid, with its maximum some 
hours after noon. Their observations were obtained in late 
1985, a time of minimum solar activity. 

In the method presented here we used GPS measure- 
ments that provide TEC values as a function of time. Such 
GPS-based TEC determinations were first successfully ap- 
phed to geodetic VLBI by Sardon et al. ( |1992| ). 



time of broadcasting of the codes P or C /A with the time 
of receipt of these codes. An observable that corresponds 
to dT can be obtained by cross correlation of the signal 
received from the satellite at each band with a reference 
signal generated at the receiver and "tied" to GPS time. In 
this case, dT — 4)/2tiv, where is the total phase change 
of the signal during propagation from satellite to receiver 
and V is the center frequency of this transmitted signal. 
This observable is the phase delay of the carrier signal. 

The GPS observable is affected by the following: the 
ionosphere, the troposphere, 27r phase ambiguities (equiv- 
alent to multiples of of 634.75 ps and 814.60 ps, respec- 
tively, for LI and L2), multipath (e.g., from signals that 
are reflected or scattered into the receiver antennas from 
objects nearby to it), different effective location of re- 
ceivers for LI and L2, instrumental delays (different for 
LI and L2), degraded coding of the signals, and clock er- 
rors. The contributions to the observables of the largest 
of these effects can be sharply reduced by application of 



suitable techniques (see, e.g., Blewitt 1990) 



2. The Global Positioning System and the TEC. 

An introduction to the Global Positioning System (GPS) 
can be found in Hofmann-Wellenhof et al. (1997). A 



main use of the GPS system is to determine the position 
{x,y,z,t) of a GPS receiver on Earth's surface. The sys- 
tem consists of a constellation of 24 satellites broadcasting 
electromagnetic signals in two narrow frequency bands, a 
set of monitoring ground sites, the Master Control Sta- 
tion, and GPS receivers. The 24 satellites orbit the Earth 
in near circular orbits with a 12 hr period, at a height of 
about 20,200km, and an inclination of 55°. The space- 
craft are in six orbital planes with four satellites nearly 
equally spaced along the orbit in each plane. At any mo- 
ment, from any point on Earth, it is possible to detect 
signals simultaneously from 7 to 9 of these satellites. Each 
satellite broadcasts a block of data every 30 seconds, con- 
sisting of a description of its orbit and of GPS time, as 
well as a pseudo-random code every millisecond (coarse- 
acquisition C/A, for civilian use) and a precision code (P, 
266 d period, for military use), usable to determine more 
accurately the position of the ground receiver. 

The oscillators of the satellites generate a fundamental 
frequency — 10.23 MHz, which is the P-code frequency. 
The C/A-code frequency is Vq/IO. The GPS signal is emit- 
ted at two frequencies, 154 j/q and 1201^0 (LI and L2, re- 
spectively, 1,575.42 MHz and 1,227.60 MHz, or AA19 and 
24.4cm). L2 carries only the P signal, and LI both P and 
C/A signals. 

2.1. The GPS observables. 

The main GPS observable is dT, the time of transit of 
the signal from the satellite to the ground receiver. The 
value of dT can be determined, in effect, by comparing the 



2.2. Obtaining the TEC from GPS data. 

The GPS observable can be modeled as a function of dis- 
tance from satellite to receiver, ionospheric delays, tropo- 
spheric delays, clock errors, and i nstrum ental phase- and 
group-delay biases (Sardon et al. 1994). The ionospheric 
term can be estimated by a combination of the LI and L2 
observables. We can denote the TEC for any observation 
direction as P^{t). This TEC is defined along the line of 
sight from the radio telescope k to the radio source i and 
can be expressed approximately as a function V (verti- 
cal value of TEC) of time t and the intersection point P"^ 
( "ionospheric point" ) of the line of sight from k to i with 
the (average) F2 layer of the ionosphere (at an altitude 
of hp^ = 350 km), times the obliquity or slant function 
5(e5j), defined as the secant of the zenith angle at the 
ionospheric point (see below for geometry clarification), 
which is a function of the elevation angle ej. of the obser- 
vation: 



lUt) = S{el)-V{Pl,t). 



(4) 



The positions of the involved sites and the ionospheric 
points can be expressed in a geocentric coordinate system 
{X,Y,Z) or in a geocentric-solar one {'^,x,Z') with the 
Z' axis directed toward the Sun from the Earth center, 
the angle in the XY-plane (measured counterclockwise 
from X) and x the angle with apex at the Earth's center, 
measured from the direction to the Sun {Z') to the direc- 
tion to Pj. This latter coordinate system is useful since 
the ionosphere is roughly time-independent in this refer- 
ence frame. Following Sardon et al. ( |1994D , we replace V 
for a GPS site j and satellites I by its locally linear ap- 
proximation in Pjj (ionospheric point towards satellite I) 
using the (^',Xi ^')"COordinate system: 

lUt) = S{e\y[A,it)+B,{t)dmt)+C,{t)dx\]+K'+K,.{5) 
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Here, the coefficients for each site are A (0'^ order, verti- 
cal), B {V^ order, ^-direction), C {V^ order, x-direction). 
We have also introduced the instrumental GPS satellite 
K'^ and receiver Kj biases. Aj{t) is the vertical TEC at 
site j, d^'' (t) = — ^j, and dXj- = ~ Xj are the co- 
ordinates 5" and X of ionospheric point of the line of 
sight from j towards the satellite I minus the correspond- 
ing coordinates of the GPS site j. The coefficients A, B, 
C, can be determined, and the if-biases largel y rem oved, 
using a Kalnian filtering method (Herring et al. 1990| ) with 
the data from different satellites / (about 8 at a time) to 
obtain an estimate of the TEC from GPS data. The K 
biases can be due to such effects as errors in the estima- 
tion of phase ambiguities and multipath (see Sardon et al. 
f994 ). In sum, we produce estimates of the total electron 



content of the ionosphere using data from a network of 
GPS stations and a multiplicity of satellites. 



Zenith angle = arcsin(H cos 
Ionospheric point / 

, w Towards / 




Site k 



Earth angle £1 
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Fig. 2. Geometry of the obliquity or slant factor, which 
is the secant of the zenith angle, defined by the elevation 
el (see text for definition of S). Curvatures of the Earth's 
surface and the ionosphere F2 layer are exaggerated for 
clarity. 



3. Vertical TEC from GPS and TEC for VLBI. 

Here we describe the formulas we used to estimate the 
TEC along the paths of a VLBI observation from the val- 
ues of the vertical TEC (evaluated by the method de- 
scribed in Sect. 2.2) at a GPS site near each one of the 



VLBI sites. 

Consider a GPS site j at latitude Q 
Xj , and a VLBI site k at position (Cfe , Afe 
The coordinates of P^, denoted by (Cp* , Ap^ 
the radio source is at elevation e 



and longitude 
) (see Fig. i. 
) at the epoch 



(Klobuchar 1987) 



J, and azimuth a^, are 



^pi = arcsin(sin cos -I- cos C/c sin £j cos a* ) 

/ sin f * sin a) \ 

Xpi = Xk + arcsin ^ , 

\ cosCp. / 



(6) 
(7) 



where {—n/2 < (k < '^f'^) ^^'^ the arcsin-function in A 
holds for values in the interval (— 7r/2, 7r/2), appropriate 
for Eq. (|) since GPS and VLBI sites are nearly collocated. 
£1 is the angle, measured from the center of the Earth 
between the line to site k and the line to the ionospheric 
point for radio source i: £1 = tt/2 — e\ — arcsin(2 cos 
(see Fig. ^, where S = Rq/{Rq + /iPa); is the Earth's 
radius, and S « 0.948 for hp^ = 350km). The local time at 
the ionospheric point is t = (Ap; /15)+UT hr (A in degrees, 
with A positive to the East). 

We assume that the vertical TEC Api (t) at the iono- 
spheric point is given in terms of Aj{t) by: 



Ap,{t) 



A, it 



15 



(8) 



This relation effects a longitude correction but ignores any 
latitude dependence of the TEC values. This approach is 
reasonable for mid-latitude stations and sources at high 
declination since GPS sites can be collocated with VLBI 
sites, or at least placed relatively near to them, and since 



TEC changes more rapidly with longitude than with lati- 
tude. 

The slant factor was defined above as the secant of 
the zenith angle, which is arcsin(S cos e).) (see Fig. |2[ and 
Klobuchar 1987). Thus, the TEC at the ionospheric point 

gives 



PI, mapped by the slant factor (Eq. ^ 

Ikit) — sec{arcsin(Scose^)} • Api{t), 



(9) 



which, with Eq. (g) yields our estimate of the ionospheric 
delays at the indicated VLBI site. The overall ionospheric 
effect on a VLBI observable is a simple linear combination 
of the effects from each of the two sites involved in the 
observation. 

4. A case to test the method: VLBI observations 
of the Draco triangle. 

Progress in high precision phase-delay difference astrom- 
etry has been made by Ros et al. ( 1999 ) through VLBI 
observations of triangle formed by the radio sources 
BL 1803-^784, QSO 1928-^738 and BL 2007-^777, in the 
Northern constellation of Draco (the Dragon). The obser- 
vations were made simultaneously at the frequencies of 2.3 
and 8.4 GHz at epoch 1991.89 with an intercontinental in- 
terferometric array. The angular separations among these 
radio sources were determined with submilliarcsecond ac- 
curacy from a weighted-least-squares analysis of the dif- 
ferenced and undifferenced phase delays. The modeling of 
these astrometric VLBI observations was sufficiently accu- 
rate to estimate reliably the "27r ambiguities" in the differ- 
enced phase delays for source separations of almost 7° on 
the sky. For such angular distances, this accurate "phase 
connection" at 8.4 GHz, yielding the phase delays with 
standard errors well within one phase cycle over the en- 
tire session of observations, was demonstrated at an epoch 
of sol ar maximu m. As in earlier works (e.g., Guirado et al. 
1995, Lara et al. 1996 ), after phase connection the effects 



of the extended structure of the radio sources were largely 
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Fig. 1. Geometry used to estimate TEC from GPS measurements (see text). The vertical lines here are radial (from 
the Earth's center) and the tangent planes drawn are approximations to the surface of spheres (Earth and F2 layer of 
ionosphere) having, e.g., the VLBI site and its ionospheric point as points of tangency. 



removed from the phase-delay observables. The effects of 
the ionosphere were also mostly removed, via the GPS- 
based method described in this paper. 

In 1991 the number of available GPS sites was small, 
and only data from Goldstone and Pinyon Flats in the 
US, and from Herstmonceux and Wettzell in Europe (see 
Table for details) were available to be used for our ex- 
periment. GPS data from the two US sites were used for 
the VLBI sites at Fort Davis (TX), Pie Town (NM), Kitt 
Peak (AZ), and Los Alamos (NM) to estimate the TEC 
for observations at these sites. Similarly, GPS data from 
both European sites were used for Effelsberg (Germany). 
The VLBI observations were carried out from 14hrUT on 
20 November to 4hrUT on 21 November 1991. In Fig. | 
we show the corresponding vertical TEC values for one 
of the GPS sites on each continent. We see the dusk and 
night part of the data for Wettzell, and the daylight data 
for Pinyon Flats. 

From the GPS-based estimates of the TEC along the 
lines of sight for each VLBI site, we calculated the iono- 
spheric contribution for each baseline and epoch of ob- 
servation. These ionospheric contributions were removed 
from the VLBI observables. For the phase delays, these 
contributions ranged in magnitude from to 1.2 ns at 
8.4 GHz for intercontinental baselines, and were less than 



Table 1. Positions of the GPS sites and the VLBI sites 



(Ros et al. 1999) and maximum and minimum values of 
the vertical TEC value on 20/21 November 1991 (see text 
for definitions of notation). 



GPS 






^max 


VLBI 




site 




[TECU] 


[TECU] 


site 




DSIO' 


35.2 N 116.9 W 


5 


57 


FD" 


30.6 N 103.9 W 


PINl^ 


33.6N 116.5W 


5 


57 




34.3 N 108.1 W 










KP= 


32.0N 111.6W 










LA"* 


35.8 N 106.2 W 


HERS^ 


50.9 N 0.3 E 


8 


38 




50.3 N 6.8 E 


WTZl'' 


49.1 N 12.8 E 


8 


39 







^ Goldstone Deep Space Tracking Station, NASA, CA, US; 
2 Pinyon Flats Observatory, SIO, UCSD, CA, US; ^ Herstmon- 
ceux, East Sussex, England; * Kotzting, IFAG, BKG, Bavaria, 
Germany. 

° Fort Davis, VLBA, NRAO, TX, US; Pie Town, VLBA, 
NRAO, NM, US; " Kitt Peak, VLBA, NRAO, AZ, US; Los 
Alamos, VLBA, NRAO, NM, US; " Effelsberg, MPIfR, North- 
Rhine- Westphalia, Germany. 



0.1 ns for continental baselines. The intercontinental base- 
line lengths range from ^--^7800 to ^8300 km, and the dif- 
ferences in local time are ^7.5 hr (see Fig. ^). Since the 
ionospheric effect is the combination of effects for both 
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Fig. 3. Values of the vertical TEC at the GPS sites of 
Pinyon Flats (CA, US) and Wettzell (Germany) (Ros et 
al. 1999). These values were used to compute the TECs 
at nearby VLBI sites. We assumed an error of 1.5TECU 
(represented by the dotted lines), much larger than the 
statistical standard errors. 



antennas, the ionospheric delay is quite important for our 
intercontinental baselines. By contrast our US-continental 
baseline lengths range from ^200 to '^750 km which cor- 
responds at most to about 30 minutes difference in local 
time and a much smaller ionospheric effect on the VLBI 
observables. 

Eq. (H) provides the usual way to largely remove the 
effect of the ionosphere on dual-band VLBI observations. 
The phase delays from our 2.3 GHz observations could not 
be freed from 27r ambiguities due to the large scatter in 
these data. However, unambiguous but less precise group 
delays were available at both 2.3 and 8.4 GHz. The iono- 
spheric contributions were estimated from these group- 
delay data. In Fig. || we show the comparison of GPS- 
based and dual-band VLBI-based ionospheric delay esti- 
mates at 8.4 GHz. We show four of the ten available base- 
lines as representative examples: two intercontinental and 
two continental-US ones. Apart from the larger dispersion 
in the group-delay than in the GPS-based data, this com- 
parison provides a good independent confirmation of the 
reliability of the GPS-based method for the correction of 
VLBI data. The error bars for the group delays shown 
in the figure are the appropriate combination of the sta- 
tistical standard errors for the data at these frequencies 
(see Eq. The statistical standard errors for the GPS 
estimates are each about 0.2TECU. We assume a much 
larger standard error -1.5TECU- to try to account for 
possible inaccuracies not estimated with the Kalman fil- 
tering, such as incorrect values for hp^j the consequent 



changes in the mapping function and in the eventual posi- 
tion of the ionospheric point, and other unmodeled effects. 
Thus, we infer a corresponding contribution of '-~^30ps to 
the standard errors of the phase delays at 8.4 GHz. For 
the data presented in Fig. ^ the root-mean-square of the 
differences between the results from the two methods are 
below 0.15ns for intercontinental baselines and 0.10ns for 
continental ones. 
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Fig. 4. Comparison between the ionospheric delays at 
8.4 GHz estimated by the combination of group delays at 
8.4 and 2.3 GHz from the VLBI observations of Ros et 



al. (1999) (points with error bars) and the corresponding 
ionospheric delays estimated by the GPS TEC determina- 
tion reported here (solid line). Data from four representa- 
tive baselines of a total of ten from the VLBI experiment 
are shown here. Root-mean-square values for the differ- 
ences are, respectively, of 0.11, 0.15, for the upper panels 
and of 0.10, 0.08 ns for the lower panels. 



5. Conclusions 

Assuming that the ionosphere can be modeled usefully as 
a thin shell surrounding the Earth at a height of 350 km, 
we estimate the TEC for the line of sight from the GPS 
site to the satellite using ground reception of GPS signals 
and knowledge of the orbital parameters of the GPS satel- 
lites. Having a GPS site near a VLBI site, we can reliably 
"transfer" this estimate to the TEC for the line of sight 
from the radio telescope to the radio source. These TEC 
values allow us to correct the VLBI observables for the 
effects of the ionosphere for any frequency. We made such 
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corrections for our phase-delay data at 8.4 GHz (Ros et 
al. 1999) for an epoch at which there was a paucity of rel- 



evant GPS data. The estimates of the ionospheric delays 
provided by the VLBI measurements of group delay at 
2.3 and 8.4 GHz differ from the corresponding delays ob- 
tained from GPS data to within root-mean-square values 
below 0.15ns for intercontinental baselines and 0.10ns for 
continental ones. Thus, we have shown, in particular, that 
the GPS determination of TEC can be successfully used 
in the astrometric analysis of VLBI observations. 

The density of the network of GPS sites has increased 
dramatically since 1991 and the accuracy of the TEC de- 
duced from GPS data has improved significantly. Given 
such progress, the approximations used in this paper are 
no longer necessary. Now GPS estimates for the vertical 
TEC are available from virtually every land location all of 
the time and thus for every VLBI observation. 

The advantages of GPS compared with geostation- 
ary beacons which used Faraday rotation to determine 
TECs are notable: global land coverage for GPS is avail- 
able from geodetic networks (e.g., the International GPS 
Geodynamics Service); the TEC estimates do not depend 
on assumptions about the Earth's magnetic field; LI and 
L2 GPS data are available over the internet in standard 
formats (e.g., RINEX: Receiver INdependent EXchange). 
We also note that the ionosphere cannot always be rep- 
resented by a thin shell model with good accuracy; more- 
over, more accurate models can be devised and suitably 
parameterized given the tomography-like sampling of the 
ionosphere provided by the GPS. The biases in the GPS 
observables have to be properly corrected in the estimates 
of TECs. Davies & Hartmann ( 1997 ) set an upper limit of 
3 TECU for the present agreement of GPS TEC with the 
results from other methods, but only of few 0.01 TECU 
for the relative errors of TEC estimates over the course 
of some hours for any given site. The latter represents 
relative changes of delay at 8.4 GHz of 0.2 ps, nearly two 
orders of magnitude smaller than the delay equivalent of a 
27r phase ambiguity. With such accuracies and the present 
network of GPS sites, the removal of most of the effect of 
the ionosphere from VLBI observations should not be dif- 
ficult, although much of the smaller-scale ionospheric ac- 
tivity cannot be adequately sampled by GPS. From dual- 
band difference VLBI astrometry and GPS data, the op- 
tical depth of the emission in radio sources can be better 
studied by comparing brightness distributions obtained in- 
dependently at each frequency band with respect to the 
same coordinate system. In sum, the introduction of GPS 
techniques should greatly improve the scientific results ob- 
tained from VLBI observations. 
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